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By weakly gauging one of the U(l) subgroups of the i?-symmetry group, Af = 4 super- Yang- 
Mills theory can be coupled to electromagnetism, thus allowing a computation of photon 
production and related phenomena in a QCD-like non-Abelian plasma at both weak and 
strong coupling. We compute photon and dilepton emission rates from finite temperature 
TV = 4 supersymmetric Yang-Mills plasma both perturbatively at weak coupling to leading 
order, and non-perturbatively at strong coupling using the AdS/CFT duality conjecture. 
Comparison of the photo-emission spectra for Af = 4 plasma at weak coupling, Af — A plasma 
at strong coupling, and QCD at weak coupling reveals several systematic trends which we 
discuss. We also evaluate the electric conductivity of Af = 4 plasma in the strong coupling- 
limit, and to leading-log order at weak coupling. Current-current spectral functions in the 
strongly coupled theory exhibit hydrodynamic peaks at small frequency, but otherwise show 
no structure which could be interpreted as well-defined thermal resonances in the high- 
temperature phase. 



INTRODUCTION 



Any thermal medium composed of electrically charged particles emits photons. The energy 
spectrum of the produced photons depends on the details of the system: the spectrum is Planckian 
when the photons are in thermal equilibrium, and deviates from it when they are not. The quark- 
gluon plasma (QGP) produced in heavy ion collisions is expected to be optically thin, because of 
its limited extent and the small value of the electromagnetic coupling a EM . Therefore, a photon, 
once emitted, should stream through the QCD plasma virtually without subsequent interaction [1]. 
In such a situation, the photon spectrum will have little to do with the black-body distribution, 
but may instead give valuable information about the properties of the medium. For instance, while 
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the experimental results [2, 3] for photon production at RHIC are currently consistent with pion 
decay plus prompt photons produced by the initial scattering of the partons from the two nuclei 
[4], there is room for photons produced in the hot plasma. 

While prompt photons really are perturbative, and pion decay photons can be calibrated from 
other hadronic signals, the most interesting signal, photon production from the medium, suffers 
from the usual problem that we only have weak coupling calculations [5] for QGP photon pro- 
duction, despite the fact that the medium is probably strongly coupled. Therefore, any guidance 
on the behavior of photon production as a function of coupling would be useful, even if it comes 
from an analogue theory which is not quite real QCD. With this in mind, we will calculate photon 
production in AA = 4 supersymmetric Yang-Mills theory, where strong coupling techniques exist. 

Consider a field theory in thermal equilibrium, and let the photon interaction with matter be 
of the form eJ^ M A^, where the electromagnetic coupling e is so small that the photons are not 
rescattered and thermalized. If T 7 denotes the number of photons emitted per unit time per unit 
volume, then to leading order in e the rate is given by [6] 

dT = " * - , rf v C<(K) , (1.1) 

7 (2tt) 3 2\k\ 1 ^ ' fco=| fc | v ; 

where 



(K) = Jd A Xe~ iKX (J™(0)J™(X)) (1.2) 



is the Wightman function of electromagnetic currents, and the expectation value is taken in the 
thermal equilibrium state. Here rj^ u = diag(— +++) is the Minkowski metric, and K = (k°, k) is 
a null 4- vector 1 whose time component is fixed by the on-shell condition k° = \k\. The Wightman 
correlator (1.2), in thermal equilibrium, is related to the spectral density, 

C<(K)=n h (k°) Xlxl/ (K), (1.3) 

where n^k ) = l/(e /3fc ° — 1) is the usual Bose-Einstein distribution function, and Xnu{K) is the 
spectral density, proportional to the imaginary part of the retarded current-current correlation 
function, 

X ^(K) = -2lmC^(K). (1.4) 

If one also adds to the theory massive fermions which carry only electric charge ("leptons"), 
then the thermal system will also emit these leptons, produced by virtual photon decay. The same 
electromagnetic current-current correlation function, evaluated for timelike momenta, gives the 
rate of lepton pair production for each such lepton species [6] : 

dTa = Wf ^°) (- K2 - 4 ™ 2 ) [-^ 2 -4-Y /2 (-^ 2 +2m 2 ) rTC<(K) . (1.5) 

Here is the electric charge of the lepton, m is lepton mass, and the correlator C< V (K) is evaluated 
at the timelike momentum of the emitted particle pair. [0(x) denotes a unit step function.] 



1 We follow the common thermal field theory convention that 4-vectors are capitalized while their components are 
lower case. 
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Expressions (1.1) and (1.5) for the production rates are true to leading order in the electromagnetic 
couplings e and eg, but are valid non-perturbatively in all other interactions. 

The electrical conductivity a of the medium is also determined by the current-current correlator, 
specifically its zero-frequency limit at vanishing three-momentum, 2 

<r = lim — rf v C<(k°, k=0) . (1.6) 

An equally valid alternative expression relates the conductivity to the the small frequency limit of 
the correlator for lightlike momenta, 

fc°->0 4J ^ l|fc|=fc° 

This form, which will be useful in our discussion of the photon production rate, follows from the 
Ward identity K^C< U (K) = combined with the diffusive nature of the hydrodynamic pole in the 
correlator. 

Both photon and dilepton production rates have previously been calculated in QCD in per- 
turbation theory to leading order in a s [i.e., up to relative corrections suppressed by powers of 
a s ] [5, 7]. There are also non-perturbative lattice estimates of the dilepton emission rate at zero 
three- momentum [8], and of the electric conductivity [9], based on attempts to fit the Euclidean 
correlator using parameterized forms of the spectral density. Reports of lattice studies of current- 
current spectral functions at non-zero momentum have appeared recently [10]. (However, the 
results of these efforts to extract real-time physics from Euclidean lattice simulations are quite 
sensitive to the assumptions made about the form of the spectral density. Assessing the reliability 
of these results is not easy; see, for example, Ref. [11].) 

In this paper, we calculate photon and dilepton production rates in SU(N C ), AA = 4 super- 
symmetric Yang-Mills (SYM) theory at finite temperature and zero chemical potential, at both 
weak and strong coupling. This theory, at non-zero temperature, mimics many features of high- 
temperature QCD. It is a non-Abelian plasma which happens to have adjoint representation 
fermions and scalars instead of fundamental representation quarks. Despite this difference in mat- 
ter field content, thermal SYM theory exhibits deconfinement, Debye screening, area-law behavior 
of spatial Wilson loops, and a finite static correlation length, just like hot QCD. But unlike QCD, 
real-time thermal properties of AA = 4 SYM theory can be studied analytically at strong coupling. 3 
Consequently, in SYM theory one may calculate, reliably, interesting physical observables at both 
weak and strong coupling. In particular, the calculation of thermal spectral functions in strongly 
coupled SYM theory is vastly simpler than the corresponding problem in strongly coupled QCD. 
Full spectral functions of the energy-momentum tensor, at strong coupling, were calculated recently 
in Refs. [13, 14], explicitly showing that strongly coupled SYM theory behaves much more like a 
liquid, rather than a weakly interacting gas of quasi-particles. At the same time, weak-coupling 
calculations of the photon emission rates in SYM theory are qualitatively similar to those in QCD, 



2 This Kubo formula for the conductivity is more commonly written in terms of the purely spatial part of the 
correlator, a = lim fc o^ |^ (fc°, fc=0). The form (1.6) is equivalent since transversality of the current-current 
correlator implies that C^ (k ,k~0) = for any non-zero frequency. 

3 N — 4 SYM theory is a conformal field theory, whose coupling is a fixed, scale-independent parameter. The 
accessibility of the strong coupling regime in SYM theory is due to gauge-string duality [12], commonly referred 
to as AdS/CFT correspondence. Though not proven rigorously, this duality has survived an impressive number of 
consistency tests and we assume its validity. 
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and detailed comparison of the results can shed light on the degree to which thermal QCD can be 
quantitatively modeled by N = 4 SYM theory. 

Our paper is organized as follows. In section II we discuss how to couple the degrees of freedom of 
N = 4 SYM theory to electromagnetism. In section III we calculate the trace of the spectral function 
X^iK) = 7]^ u Xnv{K) in strongly coupled SYM theory for arbitrary momenta. The behavior of 
the spectral functions at small frequencies is in complete agreement with the prediction of linear 
response for hydrodynamic fluctuations of the conserved charge density. We find a finite result, 
a = e 2 N 2 T/16ir, for strong coupling limit of the electric conductivity of AA = 4 SYM theory. In 
section IV we compute the current-current spectral function in weakly coupled SYM theory, for 
both timelike and lightlike momenta, for frequencies large compared to A 2 T, which is the scale 
where hydrodynamic effects become important. (Here A = g 2 N c is the 't Hooft coupling.) This 
weak-coupling analysis generalizes the corresponding calculation for QCD performed by Arnold, 
Moore and Yaffe [5]. In the final section V we compare the photon emission spectra for SYM and 
QCD at weak coupling, and SYM at strong coupling, and discuss the relevant lessons which can be 
drawn. We find that the weak-coupling behavior of SYM is quite similar to that of QCD provided 
one compares the theories at the same values of thermal masses, rather than equal values of 't 
Hooft coupling (although SYM theory has somewhat more soft photons relative to hard photons 
in comparison with QCD). The strongly coupled theory has a greater photon production rate at 
large momentum, relative to the weakly coupled theory, but less production at small momenta 
(k <C A 2 / 3 T). The production of large mass dilepton pairs is essentially identical between the 
weakly and strongly coupled theories. 



II. COUPLING Af = 4 SUPER- YANG-MILLS TO ELECTROMAGNETISM 

The field content of J\f = 4 SYM theory consists of SU (N c ) gauge bosons, plus four Weyl fermions 
if)p and six real 4 scalars 4> pq = —(fiqp, p,q = 1, - ■ ■ ,4, transforming in the adjoint representation of 
SU(N C ). The theory has an anomaly free global SU(4) i?-symmetry, under which the fermions 
transform in the 4 and the scalars in the 6. To model electromagnetic interactions, we add to the 
theory a U(l) gauge field coupled to the conserved current corresponding to a U(l) subgroup of 
the SU(4) i?-symmetry. 

We will choose the U(l) subgroup generated by t 3 = diag(i, — ^,0,0), under which two of the 
Weyl fermions have charge ±| and two complex scalars have charge ^. The associated conserved 
current is 



tEM _ l^int _ i 

" ~~ e dAv ~ 2 



- V>2 f ^2 + E <H4 + iD^Wip] • (2-1) 

p=3,4 



A summation over the SU(N C ) group index a is implied in Eq. (2.1). The covariant derivative 
acting on the scalars involves both the SU(N C ) gauge fields and the U(l) electromagnetic potential 
Afj, (with coefficient |). However, the dependence on the U(l) gauge field, reflecting quadratic 
dependence on in the scalar field part of Si nt , does not contribute to the emission rates at 
leading-order in e 2 and can be ignored. 5 So for our purposes we can treat the electromagnetic 



4 It is convenient to regard the scalar fields as components of an antisymmetric complex matrix satisfying the reality 
condition ((P Pq )1 = \e pqrs <t>rs- 

5 In the retarded current-current correlator, this term generates an 0(e 2 ) momentum independent contact term 
which does not contribute to the imaginary part of the correlator. 
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interaction as being linear in A^, with £- mt = eJ 3 A^, where J 3 is the t 3 component of the R- 
current in pure AA = 4 SYM. We further add to the theory one or more "leptons" i which are 
fermions with electric charge ee, but with no direct interactions with any SYM fields. Hence, our 
complete Lagrange density is 

C = £sym + Ant - -A - lijp + m)£ . (2.2) 

We will refer to the theory of Eq. (2.2) as SYM-EM theory. 6 For comparisons with QCD, one 
should regard £sym as modeling strongly interacting quark and gluon fields, while A^ describes 
the photon and I represents the electron and/or muon. 

The photons and leptons, once produced, are assumed to stream through the SYM medium with 
negligible further interaction, due to a small value of e 2 . Hence their emission rates, to leading- 
order in e 2 , are completely determined by the correlation function of i?-currents, ( J„(0) J 3 (x)), 
with the expectation value taken in the thermal equilibrium state of SYM theory. The evaluation 
of this correlation function can be conducted purely within SYM theory, with no further reference 
to the EM sector. Inserting the result into Eqs. (1.1) and (1.5) will yield the photon and dilepton 
differential emission rates for AA = 4 SYM. 

The choice (2.1) for the electromagnetic U(l) current is not unique. Choosing an embedding of 
U(l) within the SU(4:)r symmetry group is equivalent to choosing a specific linear combination of 
Cartan subalgebra generators. In other words, the most general choice for electromagnetic charge 
may be expressed as 

3 

QEM = J2^Q a - ( 2 - 3 ) 
a=l 

where Q 1 = \ diag(l, 1, — 1, — 1), Q 2 = \ diag(l, — 1, 1, — 1), and Q 3 = \ diag(l, — 1, — 1, 1) are a 
convenient set of Cartan generators. Our particular choice of U (1) embedding, corresponding to 
(3± = 0, = 03 = \, gives the charged fermions and charged scalars equal magnitude charges, and 
yields an SYM-EM theory which is anomaly free. 7 It also happens to make the sum of squares of 
fermion charges equal in N c = 3, M = 4 SYM and three-flavor QCD. 8 

The SU(4)r symmetry of AA = 4 SYM guarantees that the i?-symmetry current-current corre- 
lator (J^Jt) is proportional to S ab . [H ere a,b = 1 • • • 15 are 5C/(4) Lie algebra indices, and an 
orthonormal Lie algebra basis is assumed.] Consequently, if one keeps the choice of U{\) embed- 
ding completely general, as in Eq. (2.3), then the resulting electromagnetic correlator C5, depends 
on the choice of embedding merely through an overall normalization factor of 01 + /3 2 + /? 2 (which 
is 1/2 for our particular embedding). Although it would be easy to leave the choice of embed- 
ding completely arbitrary, to simplify formulas we will use our specific choice in the next two 
sections. In the final discussion we will address the question of what overall normalization of the 
electromagnetic current in SYM-EM is most appropriate when making comparisons with QCD. 



6 We would like to stress that, unlike SYM theory, SYM-EM theory does not have a known string dual description. 

7 In an arbitrary background SU(4) gauge field, the divergence of the i?-current acquires an anomalous contribution, 
d^Jl oc d abc F^F^ c . For our chosen (7(1) embedding, d 333 =2 tr(i 3 {i 3 , t 3 })=0, so our electromagnetic current (2.1) 
is anomaly-free. 

8 In Af — 4 SYM this is |(iV c 2 — 1) = 2 (counting Dirac fermion fields), while in QCD with real-world charge assign- 
ments it is 3 x (| + | + |) = 2. 
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III. PHOTON AND DILEPTON PRODUCTION RATES AT STRONG COUPLING 



The large N c , large 't Hooft coupling limit of N = 4 SYM theory in a four-dimensional Minkowski 
space at finite temperature T has a dual description in terms of the gravitational background with 
a five-dimensional asymptotically AdS metric 

ds 2 = ^2 [-f(u) dt 2 + dx 2 + dy 2 + dz 2 ] + du 2 , (3.1) 

where f(u) = 1 — u 2 , u £ [0,1], and R is the curvature radius of the AdS space. The metric 
(3.1) describes a spacetime with a horizon at u = 1 with Hawking temperature T, and a boundary 
(where one can regard the dual field theory as residing) at u = 0. 

A method for computing the retarded correlation functions of i?-currents in the dual gravita- 
tional description was formulated in Refs. [15, 16]. Subsequently, the locations of singularities of 
the retarded correlator C r ^}(K) in the complex frequency plane were found in Refs. [17, 18], and 
the spectral function at zero three-momentum was computed in Ref. [14]. Here we shall determine 
the spectral function of the i?-currents at arbitrary momenta, and relate them to the photon and 
dilepton production rates, respectively. 

At zero temperature, the correlation function has the form dictated by Lorentz and gauge 
invar iance, 

C^(K)=P,u(K)U(K 2 ), (3.2) 

where P^fT) = ri^—K^Ky/K 2 is the usual transverse projector, and K 2 =— k 2 + k 2 . For the 
corresponding spectral function one finds 9 

N 2 

X , U (K) = -2Tm C$(K) = P^(K) ^- \K 2 \ 0(-K 2 ) sgn(fc°) . (3.3) 

At non-zero temperature, rotation plus gauge invariance implies that the correlator has the form 

C$(K) = P*(K) U T (k°, k) + P^(K) U L (k°, k) , (3.4) 

where the transverse and the longitudinal projectors are defined in the standard way as Pq (K) = 0, 
PT{K) = 0, Pfj{K) = Sij - kikj/k 2 , and P^(K) = P^(K) - P^(K). Spatial indices i,j run over 
x, y, z and k = \k\. Thus the trace of the retarded two-point function is rfCV*; = 2 n T + U L and 
the trace of the spectral function is 

x ^(fc°,fc) = -UmU T {k°,k) -2ImU L (k°,k) . (3.5) 



Except for the overall coefficient, the form of the zero-temperature result (3.3) is completely determined by Lorentz, 
gauge, and scale invariance. In zero-temperature SYM, the 7?-current two-point function is protected by non- 
renormalization theorems, and thus is independent of the coupling [19]. The overall coefficient is therefore fixed 
by the one-loop spectral function evaluated in the free theory. In the electromagnetic current (2.1), the two 
Weyl fermions can be combined to form one Dirac fermion, and for the Feynman correlator one finds H F (K) = 

4^"f^~ JV °4~ 1 (1+f )[ m (-^ 2 /M 2 ) ~ 2]. Here /j, is the MS renormalization scale, a factor of \ comes from electric 
charge assignment, and the factor of (1+|) signifies that there are two charged scalars, each contributing one 
quarter as much as a Dirac fermion. The spectral function is obtained from the relation I1 F (A') = Re II(A') + 
i sign(fc°) Im U(K). At spacelike momenta, Il F (K) has no imaginary part, while at time-like momenta one has 
to choose A; — ► k°+ie which gives ln(K 2 —iek°)/fi 2 — in + ln(—K 2 )/fi 2 (for positive fc°). The spectral function 
Xvv(K) = —2P/J.V ImH(K), at large iV c , is then given by the result (3.3). 
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Both H T and T1 L contribute to the dilepton rate, but only H T contributes to the photon emission 
rate, because the longitudinal part must vanish for lightlike momenta (otherwise the correlator 
would be singular on the lightcone). According to the gauge/gravity duality prescription [12], 
two-point functions of conserved currents in SYM theory are calculated by analyzing linearized 
perturbations of a U(l) gauge field Ac (having nothing to do with the electromagnetic potential 
discussed in the previous section) on the five-dimensional AdS-Schwarzschild gravitational back- 
ground (3.1). These perturbations obey Maxwell's equations, 8a{\/ —g g AB g CD Fbd) = 0, where 
gAB is the metric of the background spacetime (3.1), and Fbd = <9bAd — 8dAb is the Maxwell 
field strength. The Bianchi identity for Fbd then implies that the electric fields Ei = F$i obey the 
equations [18]: 



w 2 f , w 2 - q 2 f 



E'L + J -E' ± + W q 2 J E± = 0, (3.6a) 



E'i + ., ^ J 2f , E[ + ?—±L E l} = , (3.6b) 
" f(w 2 - q 2 f) 11 uf 2 " 

where w = k°/{2nT), q = k/(2nT), primes denote derivatives with respect to u, and the subscript 
refers to the component which is either perpendicular or parallel to the direction of the three- 
momentum k. The equations (3.6) have singular points at it=±l, and oo. 10 At u=l (the horizon), 
the exponents are These two exponents correspond to two local solutions representing 

waves coming into or emerging from the horizon. To compute the retarded correlators, one has 
to impose the incoming wave boundary condition at the horizon, thus choosing —iw/2 as the 
correct exponent [15]. At u=0 (the boundary), the exponents for both equations (3.6) are and 
1. Solutions to Eqs. (3.6) satisfying the incoming-wave condition at the horizon can be written as 
a linear combination of two local solutions near the boundary, 

E i (u)=A i Z}(u)+B i Z}\u), (3.7) 

where the index i labels the components of the electric field (and no summation over i is implied). 
The solutions Z\ and Z] 1 are given by their standard Frobenius expansions [20] near u=0, 

Z](u) = 1 + KZf(u) In it + 6^ « + ■■-, (3.8a) 
Zf{u) = u (1 + u + fog u 2 + ■■■). (3.8b) 

All the coefficients {b^} (except b^), {b^j}, and hi are determined by the recursion relations 
obtained by substituting the above expansion into the differential equations (3.6); for example, 
hi = q 2 — w 2 . Without loss of generality, one can set bty = 0, thus fixing the definition of Z\(u). 

The correlators are essentially determined by the boundary term of the five-dimensional on-shell 
Maxwell action [15, 16, 18] 



n 2 t 2 r d A K r f f 



(3.9) 



The longitudinal equation (3.6b) also has an integrable singularity, with exponents and 2, at u 2 — 1 — w 2 /q 2 . 
When integrating the equation for spacelike Minkowski momenta, this singularity may be avoided by making an 
infinitesimal Wick rotation. 
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Applying the Lorentzian AdS/CFT prescription [15], one finds 11 

IIW k) - lim E '^ K) rfW k) - lim SJMQ (3 10) 

Choosing, for convenience, the three momentum k to lie along the z direction, so that E± = 
(E x , E y ), and E\\ = E z , and using the expansions (3.7) and (3.8), the retarded correlation functions 
reduce to [18] 

To evaluate the correlators, one isolates the incoming wave part of the fluctuating field by finding 
solutions to Eq. (3.6) of the form 12 

Ei(u) = (1 - n)-™/ 2 (1 + u)~ w I 2 yi (u) , (3.12) 

where yi(u) is regular at u = 1. Given the solution [obtained by integrating Eq. (3.6) numerically, 
if necessary], one may extract the coefficients A% and Bi from the near-boundary behavior, and 
obtain the resulting correlators from Eq. (3. 11). 13 



A. Lightlike momenta 



For light-like momenta, w=q, inserting the ansatz (3.12) into the transverse electric field equa- 
tion (3.6a) produces a hyper geometric equation, and yields the analytic solution 

E x {u) = (1 - u)- lw / 2 {l + u)- w ' 2 2 F 1 (l - \(l+i)w , -\(l+i)w ■ 1-iw; §(l-u)) , (3.13) 

where 2^1 (a, b; c; z) is Gauss's hypergeometric function. To extract the imaginary part of the 
retarded correlator, which is all we need, it is convenient to convert (3.10) to the form 

N 2 T 2 

Imn T (fe°,fc) = c — Im [f(u)F ± (u,K)*F' ± (u,K)] , (3.14) 

with F±(u,K) = E x (u,K)/E x (0,K). This formula for Im I1 T reduces to (the imaginary part of) 
Eq. (3.10) in the limit u — > 0, but expression (3.14) [which is effectively the Wronskian of F± and 
its complex conjugate] is actually independent of u. Instead of taking the u — > limit, it is more 
convenient to evaluate this in the limit u — > 1. Using the solution (3.13), we find 

ton r (J0 = _^ to *L, (3,5) 



A contact term, proportional to K 2 , is to be discarded in this expression. This contact term is real, and does not 
contribute to the physically relevant spectral function. 

The factor (l-u)- lw/2 is dictated by the incoming-wave condition at u=l. Separating a factor (1+w) w ' 2 in 
addition is a matter of technical convenience. 

Instead of integrating Eq. (3.6), with boundary condition (3.12) outward from the horizon to the boundary, and 
extracting the coefficients Ai and Bi from the near-boundary behavior, improved numerical stability may be 
obtained if one also integrates inward from the boundary to find directly the solutions Z\ (u) and Z\ T (u) with the 
prescribed boundary behavior (3.8). The coefficients Ai and Bi in Eq. (3.7) may then be determined from the 
values and derivatives of these three solutions at an arbitrary interior point within the interval [0,1]. 
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0.2 



1 2 3 4 5 

w = fc°/(27rT) 

FIG. 1: Trace of the spectral function for lightlike momenta divided by frequency, •q ,J '' J Xnv{w=q)/w, in units 
of iiV c 2 T 2 , plotted as a function of frequency, with w = k°/(2irT) and q = \k\/(2nT). At small frequency, 
X tM IJ/ (w=q) / w approaches a constant limiting value, while at large frequency X ij 'u,{ w= q)/ w fsJis as w _1 ^ 3 - 
The solid (red) line shows the exact result (3.18) while the dashed lines show the low- and high-frequency 
asymptotics (3.19). 

where the denominator is a product of two hypergeometric functions, 

D(w) = 2 Fi(l + \{l+i)w ,±(l+i)w;l+iw;fj 2 Fi(l- \{l+i)w ,-\{l+i)w;l-iw;fj . (3.16) 

With the help of the identity 2F1 (a, b; c; z) = (1— z)~ a 2F1 (a, c— b; c; z/(z— 1)) [21], the denominator 
can be written as 



D(w) =4| a f , 1 (l - + ±(l-i)w;l-iw;-l)\ 2 . 

Therefore, the spectral function for light-like momenta is 

N 2 T 2 w 



X* Jk l) =k) = -4lmU T (k°=k) 



2F1 (1 - \{l+i)w, 1 + \{l-i)w; 1-iw; -1) j 



(3.17) 



(3.18) 



This result shows that the trace of the spectral function x^(K) 1S manifestly positive, as 
it should be. (Note that for light-like momenta, xtt = Xzz, an d therefore x^n = ^Xxx-) The 
asymptotic behavior for small and large frequencies is 



[\N 2 T 2 [ w -^w 3 + 0(w 5 )], 



X f "u(w=q) 



w < 1; 



(3.19) 



k \N 2 T 2 w 2 / 3 3 5 / 6 r(|)/T(i) + O(l) , w » 1 . 



A graph of the trace of the spectral function at light-like momentum, together with the asymp- 
totics (3.19), is shown in Fig. 1. The leading small-frequency behavior agrees with that found 
earlier in Ref. [16], where it was used to evaluate the diffusion constant of i?-charge in SYM theory. 
The expression (3.18) for the spectral function is valid to leading order in the limit of large N c and 
large 't Hooft coupling. This result shows that the photon production rate for M = 4 SYM theory 
approaches a finite limit as A — ► 00. 



14 These asymptotics are derived in appendix A. A simple approximation which is asymptotically correct and accurate 
to better than 2% for all frequencies is X^(.w=q) « \N*T 2 w(l + ^ 3 5/6 r(§)/r(§)] 2 V) 1/24 /(l + \n 2 w 2 ) 1/4 . 



B. Timelike and spacelike momenta 



^X T (w,q) 



A X L (w 7 q) 





FIG. 2: Transverse and longitudinal spectral functions for time-like momenta, shown in the (w, q) plane. 
Axes are q± = w±q; the tip of the light-cone is on the left. The graphs show finite-temperature contributions 
to x T = Xxx + Xyy (left): an d X L = —Xtt + Xzz (right), plotted in units of \N 2 T 2 . The subtracted zero- 



temperature contributions are x T ( w j q)\t=o ■ 
is zero on the light-cone because Xtt(k°=k) 



tt (w 2 -q 2 ), and x (w,q)\ T =o = § {w -q 2 ). Note that x («>, q) 



Xz 



,,(k°=k). 



At time-like momenta, both I1 T and H L contribute to x^(k°,k). The mode equations (3.6) 
cannot be solved analytically for arbitrary frequency and wavevector, so we determine the spec- 
tral function numerically, as explained above. The result is shown in Fig. 2, where we plot the 
temperature-dependent portion of the transverse and longitudinal contributions to the spectral 
function in the (k°, k) plane, defined as 

A X T (k°,k) = P^(K)[ X ^(k ,k)-x^ = o( k °^)] 

N 2 

= X xx (k°, k) + X yy(k°, k)--^ [(k°) 2 -k 2 ] e((k°) 2 -k 2 ) sgn(k°) , (3.20) 



A x L (k°,k) = P^(K)[x^(k°,k)- X ^ =0 (k°,k)] 



N, 



Xz Z (k°, k) - X tt(k°, k)-j^ [(k°) 2 -k 2 ] e((k°) 2 -k 2 ) sgn(k°) 



(3.21) 



The complete result for x^ is plotted in Fig. 3 as a function of frequency for several values of the 
spatial momentum. As these plots show, x^u rapidly approaches the zero-temperature curve as the 
frequency increases. The oscillatory finite-temperature deviations are shown directly in Fig. 4. 

Slices of transverse and longitudinal spectral densities at fixed frequency, plotted as a function 
of spatial momentum, are shown in Fig. 5. Note that the longitudinal spectral density is not always 
positive (for positive frequencies). The components Xzz and xtt are individually both positive (for 
positive frequencies), but their difference xl = Xzz — Xtt can have either sign. As one moves deeper 
into the spacelike region, the spectral densities rapidly decrease. 

At small frequency and small momentum, the longitudinal spectral density xl = Xzz — Xtt has 
hydrodynamic structure which cannot be resolved in Figs. 3 and 5. The time-time and longitudinal 
space-space components should behave as 

2ujDk 2 _ , 2u?D 



Xtt(u,k) 



uj 2 + (Dk 2 ) 2 



Xzz(u,k) 



2 + {Dk 2 ) 2 



(3.22) 
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FIG. 3: Spectral function trace x^(fc°, k) (left) and x^(k°, k)/w (right), in units of N*T 2 /2, plotted as a 
function of w = k° /(2irT). The different curves correspond to differing values of the momentum; from left 
to right, q = k/(2-KT) = 0, 1.0, 1.5. The dotted black lines show the zero-temperature result. 




FIG. 4: Deviation of the spectral function trace X /i At (fc°, k) from its zero temperature limit, in units of 
N 2 T 2 /2, as a function of w = k° /(2irT). The different curves correspond to differing values of the momen- 
tum; q = k/(2nT) — (blue), 1.0 (red), 1.5 (green), and 2.0 (brown). The curves at non-zero momentum 
have cusps on the light cone [at w = 1, 1.5, and 2, respectively] where the zero-temperature result first turns 
on in a non-analytic fashion. The dotted line showing the envelope of the cusps is the plot of x^u on the 
lightcone. 



where to = k°, D is the i?-charge diffusion constant, and 3 = /3(Q 2 )/(volume) is the charge sus- 
ceptibility. As the spatial momentum k —* 0, Xtt(w, k)/co approaches a delta-function in frequency 
(times 27r5). This behavior is shown on the left in Fig. 6. The longitudinal space-space spectral 
function Xzz , divided by w 2 , is displayed on the right in Fig. 6 as a function of frequency for various 

2 

values of momentum. The Ward identity for the correlator implies that Xzz{u, k) = ^ Xtt(u, k), so 
Xtt and Xzz contain exactly the same information. But with this scaling, one sees both the diffusive 
hydrodynamic peak at small frequency for the low momentum curves, together with the approach 
of all curves to a common high frequency value of jN 2 T 2 . This constant value is precisely the 
zero temperature result for Xzz(u,k) /w 2 . Our results are consistent, as they must be, with the 
value of the diffusion constant previously found in Ref. [16], 

D = — . (3.23) 
2vrT v ; 
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0.5 1 1.5 2 0.5 1 1.5 2 

q = fc/(27rT) q = k/{2nT) 



FIG. 5: Transverse (left) and longitudinal (right) spectral functions plotted as a function of q = k/(2wT), 
for several values of the frequency, in units of N 2 T 2 /2. The different curves correspond to differing values of 
the frequency; from left to right, w = k°/(2nT) — 0.2 (blue), 0.5 (red), and 1.0 (green). The dashed black 
lines show the corresponding zero temperature result. 




0.05 0.1 0.15 0.2 1 2 3 4 5 

w = k°/{2nT) w = fc°/(27TT) 



FIG. 6: Left: Time-time spectral density divided by frequency, Xtt(k°,k)/w, plotted as a function of w = 
fc°/(27rT), for q = fc/(2vrT) = 0.1 (blue), 0.2 (red), and 0.3 (green). In the limit of vanishing spatial 
momentum, Xtt/w approaches a delta- function in frequency. Right: Longitudinal space-space spectral 
density divided by frequency squared, Xzz(k°,k)/w 2 as a function of frequency, for q = k/(2irT) = 0.2 
(blue), 0.3 (red), 0.5 (green), 1 (orange), 2 (brown), 3 (black) and 4 (pink). One sees the diffusive peak 
for small frequency and momentum, together with the approach to the zero-temperature result at higher 
frequency. 

One may also easily extract the charge susceptibility of strongly coupled SYM, 15 

S = §iV 2 T 2 . (3.24) 

Finally, in Fig. 7 we plot x% as a function of q + = (fc°+fc 3 )/(27rT) for various values of q~ = 
(k°—k s )/(27rT). The q~ = curve corresponds to light-like momenta; this curve is an odd function 
of q + . As one moves away from the lightcone by increasing q~ , the curves with q~ small compared 
to 1 clearly show hydrodynamic "wiggles" at small q + , but this structure broadens and becomes 
washed out at larger values of q~ . 



Ref. [16] found that x^{^,k) = (N^T/8tt)uj, and Xtt(w,Jfc) = (N^T/8n) ujk 2 /[u 2 + (Dk 2 ) 2 ] with D = l/2irT. 
Comparison with the form (3.22) immediately gives the stated value of the susceptibility which is, of course, 
consistent with the Kubo formula DS = lim^^o § UJ ~ 1 X fl ^{ LU J k=0). 



12 



-0 .75 



-1.5 -1 -0.5 0.5 1 1.5 

q+ = (k° + fc 3 )/(2vrT) 

FIG. 7: Spectral density trace as a function of q + = (k° + fc 3 )/(27rT) for various values of q~ = (k° — 
k 3 )/(2nT), plotted in units of N^T 2 /2. The left side of the plot, q + < corresponds to spacelike momenta, 
while the right side, q + > 0, is the timelike region. The various curves (from bottom to top, except near the 
origin) correspond to q~ = (black), 0.1 (blue), 0.2 (red), and 0.3 (green), 0.5 (orange), and 1.0 (brown). 



C. Electrical conductivity 

The electrical conductivity a can be computed by using the Kubo formula (1.6) expressing the 
conductivity in terms of the zero- frequency limit of rf' v Cp V {K'), or equivalently the zero- frequency 
slope of the trace of the spectral function of electromagnetic currents (times T). The small- 
frequency behavior of the spectral function of i?-currents in strongly coupled SYM theory was 
analyzed in Ref. [16]. Inserting the limiting low frequency behavior found in that work into the 
Kubo formula (1.6), gives 

a = e 2 ^ . (3.25) 

Using the low frequency behavior of the spectral density (3.19) for null momenta to evaluate the 
lightlike Kubo formula (1.7) yields the same value, as it must. 

The result (3.25) demonstrates that the conductivity is finite and coupling-independent in the 
limit of large coupling. Note that a is sensitive to the total number of degrees of freedom in 
the theory, and therefore is not directly useful as a means of comparing transport properties in 
different theories. A more "universal" quantity is obtained by dividing the conductivity by the 
charge susceptibility (3.24), giving 

a 1 , 
= , 3.26 

in strongly coupled SYM. This is precisely the diffusion constant of ii-charge [16], showing the 
consistency of the Einstein relation <r/(e 2 H) = D. 



D. Thermal resonances? 



In a confining theory like QCD, the spectral density of the zero-temperature current-current 
correlator will have delta-function contributions from mesons like the p and J/ip, plus narrow 
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peaks from other hadronic resonances (with widths vanishing as N c — ► oo). The delta functions 
will acquire thermal widths (which are also 1/N C suppressed) at non-zero temperature, but for some 
range of temperatures one will see easily recognizable resonances with widths small compared to 
their energies. There is some evidence from lattice studies that the J/i/j remains a well-defined 
resonance even at temperatures of a few times T c [22-25], and there has been recent discussion of 
possible signs of other "bound states" above T c [26]. 

Since it is a conformal theory, J\f = 4 SYM has no particle spectrum and the zero temperature 
current-current spectral density (3.3) is featureless. However, one may mock up a confining theory 
by considering deformations of the gravitational description of N = 4 SYM in which one cuts off the 
AdS space at some value of u = u c , so the coordinate u ranges from u = (the boundary) to u = u c 
(the cutoff). According to AdS/CFT duality, string theory on this cut-off geometry should describe 
a large N c field theory with a mass gap, and hence a discrete spectrum of bound states, determined 
by the eigenvalues of the corresponding wave equations. 16 When the field theory is considered at 
non-zero temperature, there will be a confinement /deconfinement transition at a (non-zero) critical 
temperature T c . For temperatures below T c , the relevant dual gravitational geometry remains the 
same as at zero temperature (but with time periodically identified when analytically continued 
to Euclidean signature). Above T c , the field theory will be in a deconfined plasma phase, and 
the appropriate dual geometry is AdS-Schwarzschild (times 5 5 ). In the simple hard- wall model, 
properly comparing the gravitational action of these geometries (which determines the free energy 
of the thermal field theory) shows that the transition to the AdS-Schwarzschild geometry occurs 
when the hard- wall cutoff is inside the horizon [30]. 

In this model, the spectral function of i?-currents in the low-temperature phase is temperature 
independent, and equal to a sum of discrete delta functions, whose locations are determined by 
the energies of the bound states, which are eigenvalues of normalizable fluctuations in the cut-off 
AdS geometry. Above the critical temperature, the hard-wall cutoff is hidden by the horizon, 
and is entirely irrelevant to physics outside the horizon. The spectral function of i?-currents is 
precisely the same as in pure AA = 4 SYM. The results of section IIIB explicitly show that the 
spectral functions have no structure which could be interpreted as peaks corresponding to narrow 
resonances which survive in the high-temperature phase. Thus in the hard-wall model, bound states 
"dissolve" completely at the confinement/deconfinement transition. Whether this reflects physical 
features which may be shared by real QCD, or is just a pathology of the hard- wall AdS/CFT 
model, is not completely clear. We suspect it is a generic feature of light hadrons in confining 
gauge theories at iV c = oo. 

IV. PHOTON AND DILEPTON PRODUCTION RATES AT WEAK COUPLING 

A. Dilepton production 

When A is sufficiently small, one may use weak-coupling methods to compute the photon and 
dilepton emission rates. The easiest process to analyze perturbatively is the dilepton production 
rate, because it arises already at 0(A°). The simplest way to structure the calculation is to 
compute r/^ u C< u (K) directly. This requires evaluating the single cut, one loop graph shown in 
Figure 8. The cut lines have the propagator replaced by the appropriate statistical function times 



Such "hard- wall" cut-off models were discussed from the earliest days of AdS/CFT, and represent the simplest 
version of how confinement may be realized in the dual gravity description [27-29]. 
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FIG. 8: Graph needed to compute the dilepton production rate to lowest order. The diagonal slash represents 
a "cut" through the diagram; the solid line is either a fermion or a scalar. 

the discontinuity in the propagator (the difference between +ie and — ie prescriptions), which at 
this order means the substitution of —i/(p 2 +m 2 ) by 2ir n(p°) 5{p 2 +m 2 ). Noting that the sum of 
the charge squared for all SYM Weyl fermions coincides with that for the charged scalars, and 
equals \{N 2 — 1), the leading order result for the Wightman function is 17 

r,^C<(K) = -i(iV c 2 -l) J ^2n5(P 2 )2n5((K-P) 2 ){n h (p )n h (k -p°)(2P-K) 2 

+ n i (p )n i (k°-p ) Tr [\ (l- 7 5 )f 7 ^(^-f ) 7/i ] } 
= -i(iV c 2 -l) J P 2d P d ™*°^ S (K 2 + 2pk° - 2pkcos9 pk ) 

x K 2 {n h (p°)n h (k°-p ) + 2n { (p°)n i (k°-p )} , (4.1) 
where we used (P—K) 2 = to set 2P ■ K = K 2 . The integrals are straightforward and give 



rTC<{K) = 



(-K 2 ) 
16-7T k 

(N 2 -l) 



/ o _ 2 fc dp [2m (p) n { {k°-p) + n h (p) n b (k°-p)] 



16vr 



(-K 2 )n h (k°) 



2T l+ e -( fc() - fe )/ 2T T l- e -( fe0+fe )/ 2T 
3 ~ T ln i+ e -(fc°+fc)/2r + ~k ln i_ e -(fe°-fc)/2T 



, (4-2) 



for timelike K, K 2 < 0. Inserting this result into Eq. (1.5) yields the actual dilepton emission 
rate. To express this in terms of the spectral weight, one must merely remove the factor n^(k°); 
the result differs from the vacuum result by the two logarithmic factors inside the square bracket, 
which vanish exponentially for (A; — k) large compared to T. 



B. Photon production from 2 <-> 2 scattering 

The one-loop result (4.2) for rj^ v 'C^ V {K) , which is independent of A = g 2 N c , vanishes on the light 
cone, K 2 = 0. Consequently, the spectral weight for lightlike momenta first arises at the two-loop 
level. Physically, the timelike spectral weight represents the splitting of a timelike virtual photon 
into a pair of charged particles, a process which occurs even in the absence of strong interactions. 
In contrast, the lightlike spectral weight represents real photon production, which only occurs via 
scattering processes and therefore involves powers of A and higher loop orders. A complication is 
that in a thermal system, the expansion of physical quantities in powers of A is not the same as a 
diagrammatic expansion in the number of loops. This is a consequence of sensitivity to energy and 



The weak-coupling results of this section are valid for arbitrary N c ■ 
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FIG. 9: Relation of the 2-loop contribution to the current-current correlator and tree-level 2^2 scattering 
diagrams. Interference contributions arise when the gluon runs from one solid line to the other (not shown). 

momentum scales which are parametrically small compared to T. For lightlike momentum, this 
complication arises at the first nontrivial order, and requires an infinite resummation of diagrams to 
find the leading order weak-coupling photon production rate [31]. This rate can be understood as 
the sum of a contribution from Compton-like 2^2 scattering processes [32, 33] and near-collinear 
bremsstrahlung and pair-annihilation processes [34], which are further corrected due to the LPM 
effect [35-37]. A complete treatment for the QCD plasma is given in Refs. [5, 31], and we will 
extend it here to the case at hand. The main new complications are the appearance of scalar fields 
and Yukawa couplings in the 2 <-> 2 processes, and the addition of bremsstrahlung from charged 
scalars, which fortunately was already discussed in Ref. [31]. 

Consider first the 2^2 particle processes where two SYM excitations collide to produce an 
SYM excitation and a photon. They arise in the current-current correlator at the two loop level, as 
illustrated in Figure 9. Calculation of these contributions involves integrating the squared matrix 
element for each possible production process over all possible momenta of the SYM particles, with 
appropriate population functions and Pauli blocking or Bose stimulation functions on final states. 
The resulting contribution to the photon production rate has the form 

e 2 V^C<^ 2 (K) = J **f/ 2 *M (2 7 r)^ 4 (P + P'-^-if')El^| 2 ^^K(y)[l±^(^)], 

acd 

(4.3) 

where p,p' represent the momenta of incoming particles of type a, c, k' is the momentum of an 
outgoing particle of type d, n a = rif or according to the statistics of species a, and the ± sign is 
+ if d is a boson and — if d is a fermion. (Note that [l±n^(fc')] = e k ' l T n^k') in either case.) All 
external states can be treated as massless, since thermal corrections to their dispersion relations 
are suppressed by a power of A; therefore p° = p = \p\. The sum J^acd runs over species type, color, 
spin (including the photon spin), and particle/antiparticle where appropriate. We have computed 
these summed matrix elements for the SYM theory under consideration; the result, organized by 
the spins of external states, is presented in Table I. 

Those matrix elements with \ jt or 1/u behavior lead to small-angle divergences in the photon 
production rate. The best way to see this is to choose coordinates with the z axis aligned with k. 
One shifts integration variables from p to q = k — p and uses the spatial momentum conserving 
c)-function to perform the k' integration. Introducing a dummy integration variable to via 

1 = J duj5{uj + k - p°) , (4.4) 

the integration measure in Eq. (4.3) can be reduced to [5] 

f d?pdZp'd*k' , , 4f4/ „ . , N 1 f°° , ( m *- a -«l , f°° , . f 2n d<p 

/ 7 a n ,n Tin 2tt) 4 J 4 (P+P'-K-K') = - dq (hj dp / -- , 

J (27r)92p°2p'°2A;'° K ' K ' (4tt) 3 A; J J_ q J W y J 2vr 

(4.5) 
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Process 



FF -» 7 G 



FG -» 7F 



FF->-fS 



FS -» 7F 



SG^^S 



Diagrams 




E|M| 2 = e 2 A(jV c 2 -l)x 



2[ - + - 



f s 



2 [ 3 — + 3 — + 2 



; ! -3 ^ - 3 I - 2 



TABLE I: Fully summed squared matrix elements for all processes, organized by the spin of the participants: 
F=spin-l/2 fermion, S'=spin-0 scalar, and G=spin-1 gluon. The summation over spin, species label, color, 
and particle/ antiparticle has already been conducted; for instance, the FF — ► 7$ contribution includes 
production of both neutral and charged scalars, of either charge. For the processes involving both fermion 
and scalar lines, which two of the three diagrams contribute depends on which two of the three external 
states carry electric charge. 



and in terms of these variables, u/t ~ 2kp'(l— coscf))/q 2 at small q. (In this section, q=\q\, not to be 
confused with the normalized momentum in section III.) The 1/q 2 behavior of the squared matrix 
element makes up for the two powers of q in the dq and duj integrations, leading to a logarithmically 
divergent result. Of course, the photon production rate is not actually divergent; for sufficiently 
small q 2 , the calculation of the matrix element presented so far is insufficient and requires plasma 
corrections to the internal propagator which is responsible for the 1/t (or 1/u) behavior. These 
corrections become important when q 2 ~ XT 2 , and are referred to as Hard Thermal Loop (HTL) 
corrections [38]. The correction moderates the small q behavior of the matrix element and renders 
the production rate finite, albeit with an extra logarithmic dependence on 1/A. 

The coefficient of the log is quite easy to compute, using the above behavior of u/t and the 
quoted matrix elements for various processes. We find it by extracting the small q behavior of 
Eq. (4.3) and applying it to the region XT 2 <C q 2 <C T 2 . The resulting small q contribution to 
tTC< is 

rTC< (K) ~ 32 f ^ ^_ f q ^r A > f 27T # nf(k) VMl-cos0)e^ 

^ X(N 2 -l)T 2 ni (k) ln / gmax ~r \ 
4vr \qnnn~TVx) 
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Therefore, the log-enhanced part of the photon production rate is 



X(N 2 -1)T 2 m(k) 



47T 



lnA- 1/2 + 0(l) 



(4.7) 



We will use this coefficient to normalize all other contributions to the photon production rate, 
which can be written as 



t v C<{K) 



A(iV c 2 -l)T 2 nt{k) 



4tt 



In A" 1 / 2 + Ctatik/T) + O(VA) 



(4.8) 



Some, but not all, contributions to C to t may be extracted from Eq. (4.3). This is done by 
evaluating carefully the q 2 ~ T 2 region (using the full matrix element and phase space) and the 
q 2 ~ AT 2 region (using small q 2 approximations but including HTL corrections in the matrix 
element). The small q 2 region has already been handled in the literature [32, 33], 18 and the hard 
region can be handled by numerical quadratures integration of Eq. (4.3) after using Eq. (4.5) to 
reduce it to a triple integral. A similar set of integration variables is available for the constant and 
t/s type matrix elements, see Ref. [5]. For large k/T, the coefficient behaves as ^ln(2fc/T) plus a 
constant. It is convenient to separate this asymptotic behavior so we will write this contribution to 
Ctot(k/T) as C2^2(k/T) + |ln(2A;/T). Results for this quantity are given below in subsection D. 



C. Near-collinear bremsstrahlung and pair annihilation 




Bremsstrahlung Inelastic Pair Annihilation 



FIG. 10: Basic processes behind bremsstrahlung and inelastic pair annihilation. The gluon exchange leads 
to a small angle Coulombic scattering, and the photon is produced by nearly collinear initial or final state 
radiation. 



It is not obvious that the previous analysis [32, 33], which treated the soft momentum region in ordinary QCD, can 
be applied unmodified to J\f=4 SYM. They can, however. First, note that it is only processes involving fermions 
(quarks) which give rise to the log, which arises when the quark momentum is small. Second, the result in the 
literature only depended on the form of the quark self-energy at soft momentum — the fermionic hard thermal loop 
(HTL) self-energy. This is actually the same, up to the overall coefficient, between QCD and SYM, even though 
3/4 of the SYM fermionic self-energy comes about from interactions with the scalars. To see this, recall that the 
gluon contribution to the HTL fermion self-energy arises from the loop integration (in Feynman gauge) , 

- f K YS{Q+K)YG»„{K) = ^ r W 7 M (<g+ff)7^ {K l Q y -^- ( 4 -9) 

The gauge choice is irrelevant when we take the (Q <C K) HTL piece, since this piece is gauge invariant. The 
Yukawa interactions give rise to a loop integral of form 

E(Q) = -2^S(Q+K)A(K) = -2 — i^ ^l . (4.10) 

Using 7 M 7 a 7 M = —27™, the gluonic loop contribution immediately collapses to the same result. 
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Besides the 2 <-> 2 processes just considered, photons are also produced at leading order by 
bremsstrahlung and inelastic pair annihilation, illustrated in Figure 10. These contributions arise 
because small-angle Coulombic scattering is very efficient; the rate per particle of Coulomb scat- 
tering in the thermal medium is O(XT), rather than 0(\ 2 T) as might naively be expected. And 
as usual in gauge theories, initial or final state radiation is an efficient process; photon radiation 
occurs in 0(e 2 ) of such scatterings. This leads to a photon production rate which is 0(Ae 2 T 4 ), the 
same order as the 2 <-> 2 processes just considered. 19 

Unfortunately, the calculation of photon production via these processes is a little more com- 
plicated than just evaluating the graphs of Fig. 10. The physical reason that initial and final 
state radiation is so efficient is that the wavefunction of the radiated particle, emerging at a small 
(collinear) angle, can overlap with the emitter for a long time, so the amplitude builds up co- 
herently over this large formation time. But in a medium, further scatterings may occur within 
this coherence time; photon radiation from different scattering events can be partially coherent, as 
noted by Landau over 50 years ago [35-37]. One should therefore consider emission from a charge 
carrier as it moves through the medium, making a series of small-angle scatterings. The photon 
emission vertex can appear at any point along the trajectory; in computing the probability for an 
emission, one must integrate over this time separately in the amplitude for the process and the 
conjugate of the amplitude. Hence, there is an integral over the time difference between the photon 
vertex in the amplitude and the conjugate amplitude. Because the energy of a state with a particle 
of momentum (p+k) differs from the energy of a state with a particle of momentum p and a photon 
of momentum k, there is a phase difference which grows as the time difference becomes large. One 
must correctly incorporate the effects of multiple small-angle scatterings occurring during this ex- 
tended emission process. The photon production rate is then determined by summing the resulting 
photon production from a particular charged particle of momentum p+k over all charges in the 
medium. Leaving the detailed derivation to references [5, 31], the contribution from these processes 
to the current-current correlator boils down to 

< iV c 2 -l f°° dp ( Mk+p) [l- ni (p)][p 2 + (p+k) 2 ] n h {k+p) [l+n b (p)A 

V ^,brcm+pairl ) g J ^ ^ ^ Ap 2 (p+k) 2 2p(p+k) J 

/j2 
-!^Re[2p ± -f( P± ,p,k)], (4.11) 

where the function f(p±,p, k) is the solution to the linear integral equation 

2p± = f(p±,P,k) + J d 2 q ± [f( P± ,p,k) - f(p ± +q ± , P ,k)] . (4.12) 

In the final integral (4.11), p+k is the initial state energy of the particle radiating the photon and p 
is its final state energy; when p < the process is pair annihilation (note that [1 + n(p)] = n(—p) so 
the final state blocking/stimulation function becomes an initial state population function) and when 
(p+k) < the antiparticle is the initial particle. The difference in coefficients between the fermion 
and scalar contributions in the integral (4.11) reflects their different DGLAP kernels for photon 
emission. The integral equation (4.12) accounts for the evolution of the mixed state \p+k)(p, k\ 
through the plasma, that is, for the evolution after photon emission in the amplitude but before 
photon emission in the conjugate amplitude. It has been Fourier transformed into frequency, which 
makes it easier to evaluate but harder to interpret; the 2p± comes from the dot product of the 



Analogous processes involving scalar exchange do not have the same soft enhancement, and hence are subleading. 
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photon polarization tensor with the current; the first, imaginary term accounts for the phase due 
to the energy difference, the second term accounts for scattering events. In the imaginary term, 
m 2 ^ is the dispersion correction that a large momentum (p S> \f~XT) particle receives due to the 
thermal medium, Po n _ shell — p + m^ /2p. This turns out to be identical for scalar, spinor, and 
gauge degrees of freedom in N = 4 SYM, 

tt4 = AT 2 . (4.13) 

Gelis et al. [39] derived a very compact expression for the differential cross-section for scatter- 
ing with transverse momentum exchange q± (after integrating over the longitudinal momentum 
exchange) , 

(27r) 2 dr scatt _ mf> 

d 2 Q± qi(ql+ml)> 

where = 2m 2 XD = 2AT 2 is the static Debye screening mass. 



D. Photon production results 

The integral equation (4.12) can be solved by variational methods or by Fourier transformation 
into a differential equation. The same equation appears for both scalar and fermionic contributions 
because the scalars and fermions have the same small-angle cross-section and the same dispersion 
correction (4.13). The resulting contributions, normalized to the leading-log coefficient of Eq. (4.8), 
are presented separately in Table II as the coefficients C pa i r (from the region —k < p < 0) and 
Cbrem (from p > and p < —k), in addition to the combined value 

2k 

C tot (k/T) = \ In — + C 2 ^ 2 {k/T) + C brcm (fc/T) + C pair (£;/T) . (4.15) 

Our numerical results, within the range 0.2 < k/T < 20, are reproduced quite accurately by the 
approximate forms 

C 2 ^ 2 (x) ~ 2.01 x- 1 -0.158 -0.615 e - 0A87x , (4.16) 

and 

Cbrem (x) + C pair (x) ~ 0.954 x~ 3/2 ln(2.36 + l/x) + 0.069 + 0.0289 x . (4.17) 

The fitting form for C 2 ^2 has absolute accuracy of 0.02 in this range, and the form for Cbrem + C pa i r 
has relative accuracy better than 2%. Inserting the result for Ct Q t(k/T) into the leading order form 
(4.8) for the correlator, and then multiplying by photon phase space as shown in Eq. (1.1), yields 
the actual photon emission rate. This is plotted in the next Section. 

The various photon emission contributions are compared in Figure 11, for both N = 4 SYM and 
for three-flavor QCD. Notable differences between our SYM results and the corresponding results 
for QCD [5] include the following. 

• The function C 2 ^2{k/T) grows like T/k for frequencies small compared to T, whereas in 
ordinary QCD the corresponding growth is only logarithmic in T/k. This difference arises 
from Bose enhancement of scalar annihilation into a photon and a gluon, a process not 
available in ordinary QCD. Similarly, C pa _i T (k/T) rises at very small k/T in SYM, but rapidly 
goes to zero in QCD. This reflects pair annihilation of scalars, which is doubly Bose enhanced. 
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-0.023752 
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0.081252 





303784 


-0.044237 


1.83867 


12.5 


0.073272 





375504 


-0.057057 


2.00116 


15.0 


0.068160 





446580 


-0.065852 


2.14949 


17.5 


0.064608 





516616 


-0.073915 


2.28498 


20.0 


0.062000 





585444 


-0.077076 
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TABLE II: Individual contributions plus the combined value for the non- logarithmic constant Ctot(k/T) 
appearing in the leading-order form (4.8) for the current-current correlator (for lightlike momenta). 

• At momenta of a few times T, inelastic processes are comparable in size to the 2^2 
processes in QCD, but are relatively less important in J\f = 4 SYM. This is because the 
2^2 processes arise mostly from Compton-type scattering, which has a rate proportional 
to the fermionic thermal mass, while inelastic processes arise because of Coulomb scattering, 
with a rate proportional to the gauge boson thermal mass. In SYM the fermionic and 
gauge boson (asymptotic) thermal masses are in 1:1 ratio, while in 3- flavor QCD they are 
in 4:9 ratio. In addition, inelastic processes are suppressed in SYM by the larger thermal 
mass appearing in the first term of Eq. (4.12). They receive an extra contribution due to 
bremsstrahlung from scalars, but this is subdominant for large momentum photons because 
the DGLAP kernel coupling photons to scalars in Eq. (4.11) is less efficient at producing 
large momentum photons than the fermionic DGLAP kernel. 

• At momenta of order T or less, bremsstrahlung processes completely dominate the emis- 
sion rate in QCD, while in AA = 4 SYM the Bose-enhanced scalar 2^2 processes make a 
significant contribution down to much smaller momenta. 

Despite these difference, perhaps the most important feature is how similar the result for Ctot(k/T) 
is between the two theories. As seen in Fig. 11, the minimum value of Ctot is quite similar in the 
two theories. The growth of C to t(k/T) with increasing is a bit slower in SYM, as compared with 
Nf = 3 QCD, but in both theories the asymptotic behavior 20 is proportional to {k/T) 1 / 2 . 



This asymptotic growth is slower than linear because of the effect of multiple soft scattering (or LPM suppression) 
limiting the formation time of the radiated photon. However, quite large values of k/T are required to see this 
asymptotic behavior. 
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FIG. 11: Photon emission contributions 6*2^2 (blue dashed line), CWcm (green dot-dashed), C pa i r (red 
dotted), and C\ t (solid black) as a function of k/T. On the left are the results for Af = 4 SYM, and on the 
right the corresponding contributions for three-flavor massless QCD. 



At small frequency, for both QCD and SYM, the bremsstrahlung contribution Ch Tem (k/T) 
behaves like (T/k) 3 ^ 2 (up to a log) and becomes very large. For any fixed photon frequency, this 
is the correct leading weak-coupling behavior. However, the limits of small coupling and small 
frequency do not commute. For any given non-zero gauge coupling, this (T/k) 3 ^ 2 behavior cannot 
be valid all the way down to k = because the zero frequency limit of the correlator rj^ u C^ u (K) 
is proportional to the electrical conductivity [as shown by Eq. (1.7)], and this must be finite. Our 
treatment of bremsstrahlung and pair annihilation requires that photons emitted in response to 
a soft scattering event be nearly collinear with the emitting charged particle. This is valid for 
sufficiently weak coupling at any given photon frequency, but can fail for parametrically small 
frequency. As discussed in some detail in Ref. [31], the relevant scale at which our analysis breaks 
down is k ~ X 2 T. Below this scale, the growth of Cb rem (k/T) must be cut off and the correlator 
77^" 'C<{K) must approach a finite limiting value. 



E. Electrical conductivity 

The detailed behavior of the current-current correlator for k ~ X 2 T is hard to compute. 21 The 
analogous calculation in ordinary QCD has not (yet) been performed (though there are recent 
results for vanishing k but nonzero k° [40]). In this regime, the formation time of the photon 



21 For momenta A 2 T <fc<T the photon emission rate is determined by Eqs. (4.11) and (4.12), which (after using 
rotation invariance in the transverse plane) require the solution of a 1-dimensional integral equation in p±. The 
complication at k ~ \ 2 T is that the relevant values of p± in Eq. (4.12) become 0(T), so the approximation 
p\ <C p 2 can no longer be made. The problem then requires the solution of a 2-dimensional integral equation. For 
k <C A 2 T the angular dependence becomes trivial and the problem is again reducible to a 1-dimensional integral 
equation. 
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is so long that the emitting particle should be thought of as undergoing diffusive motion, not 
quasi-ballistic relativistic motion, during the emission event. We will not analyze this regime here. 
However, we can determine the value of the electrical conductivity, to leading logarithmic accuracy, 
and hence [via Eq. (1.7)] the limiting k — > value of the Wightman correlator. 

As already discussed, the electrical conductivity is set by the diffusion coefficient of charges. The 
diffusion length is in turn inversely related to the rate at which scatterings degrade a net current. 
The complication is that it is a functional inverse, requiring the inversion of a collision operator, 
which can be done approximately using variational techniques. We leave the detailed discussion to 
the literature [41]; here we summarize the ideas and outline the differences with respect to QCD. 

Two types of scattering process are especially efficient at scattering current-carriers, making a 
logarithmically enhanced contribution to the collision operator. The first comprises processes with 
t-channel gluon exchange ("Coulombic" processes); these have an s 2 /t 2 soft-divergent cross-section, 
cut off by plasma effects; but since the initial and final state particles carry the same charge, the 
effective scattering rate is only log divergent (i.e., logarithmically sensitive to T/mv). The second 
comprises processes involving a t-channel fermion exchange ( "Compton-like" processes) , which have 
an s jt soft-divergent cross-section and which completely re-orient the direction of the charge carrier 
if the exchanged fermion is charged. The rate of Coulombic scattering can be determined directly 
from the presentation of Ref. [41], though a new complication is that one must treat separately the 
departure from equilibrium for scalars and fermions (a complication already dealt with in Ref. [41] 
in the context of shear viscosity). The Compton-like cross-section affects both scalar and fermionic 
particles; the cross-section is 16 times the one found by naively applying formulae in Ref. [41], 
since each vertex can involve a gauge boson or one of three scalar fields; however, only half of 
these scattering processes destroy electrical current; the other half, in which a neutral fermion is 
exchanged, flip the charge carrier between a scalar and a fermion. 

Applying the technique presented in Ref. [41], taking into account the differences just described, 
one finds 

a = 1.28349 or , „ ■ (4.18) 

A 2 [ln(A- 1 /2) + o(i)] v > 

We have not evaluated the O(l) constant. Note that the 1/A 2 scaling (up to a log) is exactly 
what one would find by simply cutting off the fc -3 / 2 small frequency growth of C^ rcm (k/T) at 
k ~ A 2 T, and inserting this into Eq. (4.8). The behavior of rf ,h 'C^(K) near k = should smoothly 
interpolate between the intercept of 4Ta/e 2 and the form (4.8) which is valid for k 3> A 2 T; both the 
limiting intercept, and Eq. (4.8), should provide upper bounds on the actual value of the photon 
production rate. 



V. DISCUSSION 

Converting the differential photon emission rate (1.1) into the emission rate (per unit volume) 
as a function of photon energy gives 



dT^ a EM 

~dk~~'" : /; " 



krTC<(K). (5.1) 
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At low frequencies, the Wightman function rf v 'C^ V {K) approaches a constant proportional to the 
conductivity, as shown by the Kubo formula (1.7), and hence dT^/dk is linear in k for small k, 22 

dT^ 

—— = — 7T k . small frequency (5.2) 
dk ir z 

At high frequencies, the Wightman function is Boltzmann suppressed, as shown by the relation 
(1.3) to the spectral density. Therefore, in any equilibrium plasma, the emission rate as a function of 
photon energy must rise linearly from zero, reach a maximal value, and eventually fall exponentially. 

In weakly-coupled Af = 4 SYM theory, the hydrodynamic regime in which (5.2) applies is lim- 
ited to k Ss X 2 T. The slope crT/V 2 is parametrically large, as shown by the result (4.18) for the 
conductivity, and the maximal value of dT^/dk (which we have not evaluated quantitatively), will 
be of order aem^c — 1)T 3 (up to a log of A). For photon momenta large compared to A 2 T, the 
analysis of section IV applies and the photon emission spectrum may be expressed as 

^El = A kni{k) m 2^ [i n ( T/moo ) + C tot (k/T)] , [large frequency] (5.3) 

with the coefficient A = \(N 2 - 1) for M = 4 SYM (with our chosen charge assignments of ±?j). 
Here as before, m 2 ^ = XT 2 describes the thermal correction to hard fermion propagation in the 
medium. Note that if one ignores the k 3> X 2 T condition on the domain of validity of Eq. (5.3) 
and uses this result all the way down to k = 0, then the (fc/T)~ 3 / 2 behavior of C^ rem (k/T) will 
cause dT^/dk to be singular at k = 0, but because of the explicit factor of k in the formula (5.3) 
the singularity is integrable (and the energy-weighted spectrum is completely finite.) 

In strongly coupled N = 4 SYM theory, the photon emission spectrum is obtained by inserting 
the A=oo spectral density (3.18) into Eqs. (1.3) and (5.1), giving 



dT 7 = a EM N 2 T 3 (k/T) 2 
dk 16vr 2 e k / T -l 



( (l+i)k (l-i)k ik 



2vrT' 



(5.4) 



This is an exact expression (in the large N c , large A limit), valid for for all photon energies, 
both large and small. Equation (5.4) naturally reproduces the small-momentum form (5.2) in the 
hydrodynamic limit, with the electrical conductivity given by (3.25). Thus the slope of dT-y/dk 
at small momentum is coupling-independent in the limit of large coupling, and is parametrically 
smaller than the corresponding slope in the weakly coupled theory. The maximum of dT^/dk is 
attained at fc max ~ 1.48479 T, with the maximal rate of 

TT~1 \ 

~ 0.01567 a EM N 2 T 3 . (5.5) 



dk 



max 



At arbitrary values of the coupling, one must have (dT<y/dk) maiX = f(X)a EM N 2 T 3 , where /(A) 
interpolates between the strong-coupling result /(A— >oo) « 0.01567, and the weak-coupling max- 
imal intensity. To date, the weak-coupling expression for /(A) has not been calculated, in any 



In strongly coupled TV = 4 SYM theory, small k means fc<CT. In weakly coupled AT = 4 SYM theory, small k means 
fc<CA 2 T, which is the inverse mean-free path for large-angle scattering. One should also keep in mind that even 
though ff v C^ v {k) can in principle be computed for arbitrarily small k, the right-hand side of Eq. (5.1) ceases to 
have the interpretation of the photon production rate if k/T becomes comparable to either eN c , or (e 2 N^ /A) 2 ' 3 . 
The first constraint reflects the fact that, due to electromagnetic corrections to the photon dispersion relation, 
photons no longer propagate through the plasma like nearly lightlike excitations if k < e7Y c T. The second 
constraint reflects the scale where electromagnetic photon dispersion corrections can no longer be neglected in the 
integral equation (4.12). See Ref. [31] for details. 
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gauge theory. At large momenta, k^$>T, the photon rate in strongly coupled SYM theory decays 
as fe 5//3 e _fc//T , a stronger power than the k?/ 2 e~ k / T rate one finds in the extreme large k limit of 
the weak-coupling calculation (in which C ps .i r (k/T) ~ k 1 ! 2 due to LPM suppression). 
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FIG. 12: Photo-emission spectrum dT 7 /dk, divided by aEM^c — 1)T" 3 , in J\f = 4 supersymmetric Yang-Mills 
theory for A = oo (solid black curve), A = 0.5 (dashed blue), and A = 0.2 (dotted red). As explained in 
footnote 23, the weak-coupling curves interpolate between the rising small-frequency result (5.2), valid for 
k < A 2 T, and the falling large-frequency result (5.3), valid for k ^> A 2 T; the precise height of the sharp 
narrow peak is not known. 

Fig. 12 illustrates how the photo-emission spectrum evolves as A increases. 23 The slope at 
k = (proportional to the conductivity) decreases, the position and width of the hydrodynamic 
peak (both proportional to A 2 T) increase, and the amplitude of the spectrum for k/T J> 1 increases 
[due to the factor of m 2 ^ oc A in Eq. (5.3)]. Figure 12 also shows the strong coupling result (5.4). 
The strong coupling curve differs primarily in having the temperature T, and not some smaller 
scale, set the width of the hydrodynamic regime (in which dT^/dk is approximately linear in k). As 
a result, there is a broad maximum in the strong coupling spectrum at fc max ~ 1.5T. At sufficiently 
small frequencies, the photon production rate is largest in the most weakly coupled theory. This 
is because, at these frequencies, the photon wavelength is larger than the free path of the particles 
involved, so the charges are effectively diffusing; weak coupling means faster diffusion and therefore 
more current on such long scales. The cross-over point, below which the weak-coupling rate exceeds 
the (A independent) strong-coupling rate, scales as A 2//3 T. At large frequency, the production rate 
is greatest in the strongly coupled theory. This is because the spectrum, for weak coupling, is 
proportional to the gauge coupling A <C 1 [appearing in Eq. (5.3) as the factor m 2 ^ = XT 2 ], while 
the spectrum in the strong coupling limit has no such suppression. The results are clearly consistent 



The weak-coupling curves in Figs. 12 were generated using a smooth interpolation between the small frequency 
form (5.2) and the form (5.3) for O(l) values of k/T, with the unknown O(l) constant in the conductivity (4.18) 
set to i In | . Ref. [42] , which evaluated the complete leading-order flavor diffusion constant (or equivalently, the 
conductivity) in various QED and QCD-like theories, found that the correct constant to be added to the log equals 
this value to within ±8% for a variety of non-Abelian theories with different matter content. So this is our best 
guess for the appropriate value for SYM theory. In addition, the photon rate in strongly coupled SYM theory was 
plotted with replaced with N^—l in the expression (5.4). 
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with an expectation of smooth evolution between the weak and strong coupling regimes. 

It is instructive to compare our weak-coupling result for dT^/dk with the corresponding result for 
QCD [5]. But before this can be done in a meaningful fashion, we must first address the question of 
what normalization of the U(l) current in SYM-EM will best mimic the electromagnetic physics of a 
real QGP. One can consider various different criteria for fixing the charge (or current) normalization 
in order to compare results between different theories. 25 But for our purposes (involving the 
comparison of electromagnetic emission phenomena), the most natural criterion for fixing the 
normalization of the EM current in SYM-EM is to require that the dilepton emission spectra agree 
at large invariant mass. This amounts to demanding that the leading behavior of the current- 
current spectral density at large time-like momentum coincide in the two theories. This provides 
a simple criterion which will fix the normalization of the current, independent of the interaction 
strength in SYM, because all medium dependent corrections to the spectral density x^(-ftT) vanish 
exponentially when \K 2 \ 3> T 2 . Coinciding behavior of current-current spectral density at large 
timelike momenta also implies coinciding behavior at large spacelike momenta, and hence this 
criterion is the same as demanding that the leading short distance behavior of the ( J™(x) J™(0)) 
correlator agree between theories. (Hence this criterion is the same as that used in Ref. [29].) 

The particular choice of current we made in Eq. (2.1) happens to make the fermionic contribution 
to the large momentum behavior of the current-current spectral density the same as in QCD with 
three flavors (and N c = 3). This merely reflects the fact that the sum of squares of the fermion 
charges coincide. But the scalars of N = 4 SYM contribute half as much as the fermions to the high 
momentum spectral density. Consequently, to satisfy the condition of coinciding large momentum 
behavior of the spectral density (for Nf = N c = 3), the U (1) current used in our SYM calculations 
should be rescaled by a factor of y/2/3. This we do in the following comparisons. 

Once the normalization of the U(l) current is fixed, the only other adjustable parameter in the 
SYM emission spectra is the value of the gauge (or 't Hooft) coupling A. If one compares photon 
emission in three-flavor QCD and N = 4 SYM, at N c = 3 and the same value of the gauge coupling 
in both theories, then the SYM photo-emission rate is dramatically larger than the QCD rate, as 
shown in Fig. 13. This difference arises because the scattering rate in the SYM plasma is much 
higher than in QCD, owing to the larger number of matter fields and the fact that they are in 
the adjoint representation. In particular, the rate at which a quark undergoes photon producing 
Compton-type scattering is proportional to m 2 ^, which is 9 times larger in SYM than in QCD (for 
equal values of the gauge coupling). The rate of Coulomb scattering, important in bremsstrahlung, 
is proportional to m^, which is 4 times larger. 

Because of this difference in scattering rates, it is more appropriate to compare the two theories 
with the SYM gauge coupling adjusted to give either the same value of the Debye mass tod (in 
which case 

A SYM = 1 A QC D)) or the 

same value of for the fermions as in QCD (in which case 
A SYM _ l A QCD^ B th of these comparisons are shown in Fig. 14. With coinciding values of the 
Debye mass, shown on the left of Fig. 14, the two spectra are nearly identical at high momenta, 



The large momentum behavior of the spectral density r]^Xv(K), for lightlike momenta, differs quantitatively 
between weak and strong coupling, growing proportional to k 1 ^ 2 for weak coupling (due to C pa .i T (k/T)), and 
proportional to k 2 ^ 3 for strong coupling, as seen in Eq. (3.19). However, one can easily imagine that the true 
asymptotic behavior is a coupling-dependent power-law with an exponent u(X) which smoothly interpolates 

between the two limiting values. 

For example, requiring equality of the sum of squares of electric charges of all charged fields might seem natural. 
Or one could require equality of the EM charge susceptibility, which measures mean square fluctuations in electric 
charge density. These criteria are different (and both differ from our chosen condition). 
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FIG. 13: Comparison of leading-order photo-emission spectra for three-flavor massless QCD and Af — 4 SYM 
(for N c = 3), when the gauge coupling has the same value, a a = 0.1, in both theories. The upper solid 
curve shows the SYM result, while the lower dashed curve is the QCD result. 
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FIG. 14: Comparison of leading-order photo-emission spectra for three-flavor massless QCD (dashed curve) 
and Af = 4 SYM (solid curve) at N c = 3. Left: equal values of the Debye mass mo in both theories, 
corresponding to a SYM = 0.025 in SYM and a s = 0.1 in QCD. Right: equal values of the asymptotic fermion 
mass moo, corresponding to oj SY m = 0.011 in SYM and a s = 0.1 in QCD. 



k > 3T, while the SYM rate is larger at lower momenta (about 75% larger at k/T = 1). When the 
two theories are compared at coinciding values of the asymptotic fermion mass moo, as shown on 
the right hand plot of Fig. 14, the resulting curves are remarkably similar, with the SYM spectrum 
just a bit below the QCD result. 

A major motivation of this work was the hope that one would be able to translate knowledge 
of the photon production rate in strongly coupled SYM into a useful prediction for the rate in 
strongly coupled QCD. This seems reasonably plausible given that (for comparable numbers of 
charged quarks) the weak coupling photon spectra of the two theories are quite similar, as shown 
in Fig. 14 — provided one scales the 't Hooft coupling of SYM relative to that of QCD by a 
factor somewhere in the range 4-9. However, it is not clear if the quark-gluon plasma produced in 
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SYM, A=0 SYM, A=oo 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

\J 'w 2 — q 2 y/w 2 — q 2 

FIG. 15: Relative size of the thermal corrections to the trace of the spectral function, Ax^ik , k) divided 
by X^(k°, k)\ T =o, as a function of \K\/(2%T) = ^Jw 2 - q 2 for q = k/(2nT) = (black), 1 (blue), 2 (green), 
and 3 (red). The left panel shows the weak-coupling A — > result, while the right panel shows the strong 
coupling A — * oo result. 



heavy ion collisions can be well modeled by SYM plasma in the asymptotically strongly coupled 
A — > oo limit. A comparison of Eq. (3.18) and Eq. (4.8) at, say, k = 5T, shows that the weak 
coupling result approaches within 20% of the strong coupling result at A SYM « 4. Taking this as a 
guess for the beginning of the region where the asymptotic strong coupling result provides a decent 
approximation, and applying the above rescaling between SYM and QCD gauge couplings yields 
evP CD in the range 0.4-1.0, somewhat larger than the values of a s ~ 0.3-0.5 commonly thought to 
be relevant in the QGP formed in heavy ion collisions. Therefore, it may be best to view the photon 
production rate in infinitely strongly-coupled SYM, for k/T > 1, as an upper bound on what we 
expect the photon production rate from real QGP to be. Nevertheless, it will be interesting to see 
if incorporation of the strong coupling SYM spectral functions into models of photon production 
in heavy ion collisions improves the comparison with data. We have recently learned that efforts 
to do so are underway [43]. 

Regarding dilepton emission, our results have a simple and more positive implication. Suffi- 
ciently deep in the timelike region, thermal corrections to the spectral function become very small 
for both weak and strong coupling. This is shown explicitly in Fig. 15, which plots the relative 
correction to the zero temperature spectral density clS Si function of y/^K*. If \K 2 \ > (2ttT) 2 , 
then thermal corrections to spectral function at weak coupling are under 2%. At strong coupling 
the corrections are larger, but nevertheless no more than 15% in this regime. Since the zero tem- 
perature spectral density is independent of coupling, as discussed in Section III, this means that 
the dilepton spectrum is nearly identical at weak and strong coupling, as long as the invariant 
mass of the pair is above 2nT. This is also consistent with the modest size of the next-to-leading 
order weak-coupling result in QCD [44], which is a relative correction of — |(a s /27r)(2-7rT) 2 /i ; C 2 . 
Therefore, for large invariant mass dilepton pairs, it is undoubtedly an excellent approximation to 
use the lowest-order production rate calculation even when the coupling is strong. 

Our results for photon and dilepton emission rates at strong coupling can be extended in a 
number of obvious ways. One question is how the emission spectra at strong coupling change if 
the field theory deviates from the conformally symmetric N = 4 SYM case. As discussed above, 
perturbative rates of photon production in conformal N = 4 SYM and in (non-conformal) QCD do 
not differ significantly, when compared at the same value of the thermal fermion mass. This suggests 
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that the presence or absence of conformal symmetry does not play a decisive role in determining 
the perturbative photon spectrum. AdS/CFT can provide a similar comparison at strong coupling. 
For example, it should be possible to compute photon and dilepton emission spectra in both mass- 
deformed AT = 4 SYM [45], and in AT = 4 SYM at non-zero chemical potential, at strong coupling. 
A further question is related to the coupling constant dependence of the emission rates. It would 
be interesting to see how the spectrum shown in Fig. 12 evolves when 0(A -3 / 2 ) corrections [46] 
are taken into account. It would also be interesting to extend our analysis to theories with matter 
fields in the fundamental representation of the gauge group. Adding fundamental representation 
matter to N = 4 SYM corresponds, in the gravity dual, to the addition of D7 branes embedded in 
the AdS§ x S 5 geometry [47]. In such theories, it is natural to regard a flavor symmetry current 
of the fundamental matter fields as the electromagnetic current. Analyzing vector fluctuations 
on the D7 brane would then allow one to compute, at strong coupling, photon production from 
fundamental representation quarks of arbitrary mass added to the Af = 4 SYM plasma. 
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APPENDIX A: ASYMPTOTICS OF THE SPECTRAL FUNCTION 

For null momenta (w=q), Eq. (3.6a) for the transverse component of the electric field has the 
form 

El .^ e± + }p E± = . (A1) 

The solution obeying the incoming wave boundary condition at the horizon (u = 1) is 

E ± (u) = (l-uy^^l+uy^^F^-^l+^w^-^l+^l-iw-^il-u)) . (A2) 

We are interested in asymptotics of the solution (A2) for large and small values of w. To the best of 
our knowledge, appropriate asymptotic expansions of the hypergeometric function are unavailable 
in the literature. However, such expansions can be readily derived from the differential equation 
(Al) following the approach of Ref. [48]. 

For w S> 1, we use the Langer-Olver method [49] for constructing uniform asymptotic expansions 
(a version of the WKB approximation). Introducing new variables, 

E±(u) = J— y(u) , x = -u, (A3) 

v-/(«) 



one can rewrite Eq. (Al) as 

w 2 x — 1 

(T 



= ■ ( A4 ) 
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For w — > oo, the dominant term on the right-hand side of Eq. (A4) has a simple zero at x = and 
thus according to Ref. [49] the asymptotics can be expressed in terms of Airy functions. Moreover, 
since the coefficients of Eq. (A4) satisfy the conditions of Theorem 3.1 of Chapter XI in Ref. [49], 
one is guaranteed to have a uniform asymptotic expansion for all u G [0,1]. The asymptotic 
expansion is 

E ± (u) ~ -£SjL f- 1/4 (-u) Ai (w^a-u)) +-.., (A5) 

v-/(«) v ' 

where Ai(z) is the Airy function, 26 

C0«0 = ^3 ( i7r - 2 arCtan + lQ g ' = (1 - a 2)2 C ( a ) ' ( A6 ) 

and the ellipses denote corrections that can be systematically computed [49]. The normalization 
constant, 

C(w) = 2^e in/4 w 1/6 2~ iw - w ' 2 e i7Tw/i , (A7) 

is chosen in such a way that the asymptotic expansion (A5) coincides with the exact solution (A2) 
as u — ► 1, where E± — ► (1— u)~ lw / 2 2~ w l 2 . Using the asymptotic solution (A5), for the retarded 
correlators we find 

ttT/ , N 2 C T 2 (-t i ;) 2 / 3 3 1 / 3 r(2/3) 

n T (^=g) C — K - > y ' > , (A8) 

Correspondingly, for the trace of the spectral function we obtain 

. N 2 T 2 w 2 / 3 3 5 / 6 T(2/3) 
^=9)--^ r(1 ^ 3 - , (A9) 

In the low-frequency limit, one can solve Eq. (Al) perturbatively using w <C 1 as a small 
parameter. Since this procedure is well known (see, for example, Refs. [16, 48]), we omit the 
details. The retarded correlator for w <C 1 is given by 

n T (w=q) = -^N 2 T 2 w [l + iw log 2 - ±tt 2 w 2 + 0{w 3 )] , (AlO) 

and the resulting trace of the spectral function is 

X^=q) = \N 2 C T 2 w[l- ±n 2 w 2 + 0(w 4 )] . (All) 
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